Limit-Cycles and Chaos in the Current Through a Quantum Dot 
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We investigate non-linear magneto-transport through a single level quantum dot coupled to 
ferromagnetic leads, where the electron spin is coupled to a large, external (pseudo-)spin via an 
anisotropic exchange interaction. We find regimes where the average current through the dot dis- 
plays self-sustained oscillations that reflect the limit-cycles and chaos and map the dependence of 
this behaviour on magnetic field strength and the tunnel coupling to the external leads. 



I. INTRODUCTION 

Single-electron transport through nanostructures has 
developed into a powerful spectroscopic tool for probing 
correlations, quantum coherence, and interactions with 
the environment on a microscopic levelii^— . Some re- 
cent examples include experiments with semiconductor 
quantum dots that have provided detailed insight into 
level structures^i^. Coulomb and spin-blockade effects^, 
phonon emission^, or the statistics of individual electron 
tunnel events^. 

In this paper, we propose the time-dependent, av- 
erage current of electrons through a single level quan- 
tum dot as probe for classical non-linear dynamics and 
chaosiS. Specifically, we consider electronic magneto- 
transport through a quantum dot containing two spin- 
split levels with an anisotropic coupling between the elec- 
tron spin and an external, classical magnetic moment 
or pseudo-spin. In order to have a spin-polarized cur- 
rent through the quantum dot, we consider ferromagnetic 
leads (see Ref. [11] and references there in). 

Previous works have analyzed the anisotropic inter- 
action between two spins in a closed system under an 
external magnetic field^-"A&, showing either regular (in- 
tegrable) or chaotic (nonintegrable) classical orbits. The 
results presented here demonstrate that the signatures of 
non-linear dynamics and classical chaos of the closed sys- 
tem also persist in the non-equilibrium regime, where the 
additional coupling to the electronic reservoirs leads to an 
even richer dynamics that can be probed, e.g., by vary- 
ing the magnetic field and the tunnel rates. In particular, 
one finds a transition from a regime with damped current 
transients and a constant current, to a situation where 
the current displays self-sustained regular limit-cycle os- 
cillations, or chaotic behaviour. Limit-cycles in transport 
have also been found recently in theoretical calculations 
in mesoscopic systems coupled to mechanical degrees of 
£reedon>i^"— . 

Experimental inspiration for our model comes from the 
hyperfine interaction in quantum dots. The interaction of 
electron spins in quantum dots with surrounding nuclear 
spins is usually viewed as simply giving rise to spin re- 
laxation and decoherence2£i2i. Recently, however, trans- 
port experiments through semiconductor double quan- 



tum dots have shown non-linear current behaviour which 
has been attributed to hyperfine interaction inducing a 
dynamical nuclear spin polarization^^'- ^. The feedback 
between electron and nuclei spin polarization gives rise to 
nontrivial features in the current, including self-sustained 
oscillations^^. In this setting, the large spin of our model 
represents an effective description of the collective nu- 
clear spin system^! and the electronic part provides a 
minimal model for investigating the effects on transport 
of coupled spin-spin dynamics. 

A further potential realization of the large spin in our 
model is a magnetic impurity in a quantum dot. Several 
recent works have considered the influence of such an im- 
purity on the transport properties through the doiS^—. 
In this context, our model can be viewed as the large- 
spin counterpart of the previously studied models and 
in particular the spin-1/2 impurity model of Refs. [2^ 
and [111. This possibility is also closely related to trans- 
port through single molecular magnets^^Ti^, for which 
our large spin would map to a magnetic atom and the 
isolated levels of our quantum dot to molecular orbitals. 

We mention that our study of classical chaos in a quan- 
tum dot with coupling to an external pseudo-spin is also 
complementary to previous studies of intrinsic quantum 
chaos of, e.g., ballistic quantum dots. Those latter sys- 
tems are often analyzed with statistical tools such as ran- 
dom matrix theory^i^i. 

The outline of this paper is as follows. In Section II, 
we introduce the model Hamiltonian and the equations of 
motion. Section III presents results and a classification 
of various non-linear regimes in the form of a map in 
parameter space, and we conclude with a brief discussion 
of the experimental relevance of our finding in Section 
IV. 



II. MODEL 

A. Hamiltonian 

We investigate a quantum dot (QD) with a single or- 
bital level, coupled to an emitter (left electron lead), 
a collector (right electron lead) and to a large spin J 
(Fig. [Ta|) . An external magnetic field Bz is applied in z- 



direction which sphts the QD spin levels (Fig. Ilbp . The 
Hamiltonian for this model is 

H = Hp A + Hj + V. (1) 

Here, HpA is the Fano- Anderson model for the QD cou- 
pled to the leads, which is exactly solvable; Hj is the 
Hamiltonian for the free motion of the large spin due the 
external magnetic field; and V is the coupling between a 
dot electron and the large spin. These individual Hamil- 
tonians read: 

Hp A = ^ <^ddlda + B^_S^ 

a 

+ eikaclka^lka + {likc\i,^d„ + h.c.) (2a) 

Ika Ika 

Hj = bJ,, (2b) 
^^^^^^^ (2c) 

i—x.y.z 

where Cd is the energy of the QD level, dl/da cre- 
ates/annihilates a spin-cr electron in the dot, Si is the 
i-th component of the electron spin operator in second 
quantization, Ji is the z-th component of large spin op- 
erator, and \i is the coupling between the i-th compo- 
nents of the electron and the large spin, (^ik^jcika cre- 
ates/annihilates an electron with momentum fc and spin 
cr in lead I G {L, i?}, and 7;fc is the coupling between the 
QD and the /-th lead. Coulomb interaction in the QD 
is neglected, and thus double occupation is allowed. The 
flip-flop processes due to the spin-spin interaction are the 
origin of the non-trivial dynamics that will be shown in 
the next section. Much of the interesting dynamics found 
occurs at low magnetic fields, in particular, in a regime 
where the coupling between the electron and the large 
spin dominates Zeeman splittings {Bz ^ A). Thus, in 
this regime we believe that different ^-factors will not be 
qualitatively important, meaning the energy mismatch 
between the Zeeman splittings will not lead to suppres- 
sion of the flip-flop processes. Therefore, for simplicity, 
we assume identical g-factors for the electron spin and 
the large spin, and absorb the Bohr magneton and the 
factors into the definition of B^. 

The classical counterpart of the closed system (7/^ = 0) 
is, for zero external magnetic field (Bz = 0), a completely 
integrable system for arbitrary Aj^**, while the isotropic 
model (Aa; — \y = \z) is also completely integrable for 
finite external magnetic fields. However, in presence of a 
finite magnetic field, an anisotropic coupling between the 
electron spin and the large spin, makes the model non- 
integrable and can lead to a chaotic spin dynamics^. 
Therefore, in this work we take the coupling between the 
electron spin and the large spin to be anisotropic, and 
for simplicity we will focus on the choice 

A, - A, - A, Xy - 0. (3) 

Finally, the spin-dependent rates of the contact barriers 
are chosen so that only spin-up electrons can tunnel out 
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FIG. 1. (Colour online.) Scheme and setup of the investigated 
system, a) An electron spin S (blue arrow) in a QD is cou- 
pled via the exchange interaction A with a large spin J (red 
arrow). The QD is attached to ferromagnetic electron reser- 
voirs (brown regions), allowing electrons to tunnel through 
the QD. The large spin is isolated, b) The spin-dependent 
rates of the contact barriers are chosen so that a spin down 
electron is always trapped in the QD, while spin up electrons 
can tunnel through it (see details in the text). The large spin 
interacts with the spin of the electron trapped in the QD, al- 
lowing its spin to flip and, hence, escape form the QD into 
the right lead. 

of the QD (Fig.[lb|, and > fc^T where T is the tem- 
perature of the leads and ks is Boltzmann's constant. In 
this regime, current can flow only through the spin-up 
level of the QD. When an electron enters the spin-down 
level, it remains trapped until a spin-flip process (due 
to the interaction with the large spin) produces a transi- 
tion from the spin-down to the spin-up level, allowing the 
trapped electron to escape the QD. Notice that because 
we have taken identical g-factors for both the electron 
spin and the large spin, the spin-flip transition from the 
QD spin-down to the spin-up level conserves energy and 
the energy that the electron absorbs in the spin-flip is 
emitted by the large spin. 

B. Equations of Motion 

The equation of motion (EOM) for the expectation 
value of an operator A is 

Using this formula, we derive the EOM of each operator 
in Eq. (ga]), Eq. ^ and Eq. 
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We first observe that the length of the large spin j = | J| 
is a conserved quantity since [3^,H] = 0. Next, due to 
the interaction V, the EOMs do not close and lead to an 
infinite hierarchy of equations that needs to be truncated. 
In order to do so, we use a factorization approximation 
by invoking a mean-field approximation for V — > Vmf, 

VmF= ^^{S^{J.)+MS^)- {S^){J^)) (5) 



which neglects the term SSiSJi with 5Si = Si — (Si) and 
dJi = Ji — (Ji), i.e. the quantum fluctuations of the 
electron spin and the external spin. We expect this to be 
a good approximation when j 3> 1 and the external spin 
J can essentially be treated as a classical object due to its 
interaction with other environmental degrees of freedom. 
Furthermore, as in the semiclassical approximation we 
neglect quantum fluctuations of the large spin, we have 
no spin decay, meaning the large spin is a constant of 
motion. 

We furthermore neglect terms proportional to ^ik^i, 
namely, second order transitions due to the coupling of 
the large spin with the contacts. This is a good approxi- 
mation in the inflnite bias regime. For the electron leads, 
we perform the usual Born-Markov and flat band approx- 
imations and consider them to be in thermal equilibrium. 
Moreover, we consider the infinite bias regime, namely 
/Ui — >■ oo and ^ir — > — oo, respectively (see Appendix A 
and B for details). The resulting EOMs read 



— {n„) = X{Jx){Sy) {Sat - Sai) - T{h^) + Flo- 

I {s.) = (a(j.) + B.) {Sy) r(5,) 

^ (Sy) = -X{X){S.) + [X{l)+B,) {S.)~T{Sy) 

I {s,) = x{X){Sy) - r{s.) + i (Fit - ^Li) 

|(j.)=~(A(5.)-fi?.) {Jy) 



- {Jy) = -x{s^){j,) + (x{s,) + bA (4) 



di 

J^{l)^X{S,){jy), 



(6) 



where n„ — d^d^, F^ = F^^ + ^Ra with a =t,i, and 
F^CT and F^g. are the tunnelling rates through the left 
and right contact barriers, respectively. We have taken 
F^ = F^ = F for simplicity. In order to have current only 
through the spin-up level we take Tr]^ = 0. Therefore, 
spin- up electrons are allowed to tunnel through the QD, 
whereas spin-down electrons become trapped in it. 

The EOM for the total number of electrons in the QD 
[N — + ni) is independent of both the electron and 
the large spin components and is exactly solvable (see 
Appendix A). Thus, as 25*2 = n-f — nj,, the level occupa- 
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FIG. 2. a) Parameter space with three regions describing the 
behaviour of solutions of EOMs (Eq. (6)). Boundaries be- 
tween the regions are obtained analytically from Eq. (14). 
Region I: damped oscillations; region II: both damped and 
self-sustained oscillations; and region III: self-sustained os- 
cillations only, b) Numerically-obtained small-F region in 
the mixed region II. In the dark-coloured region, damped 
oscillations are obtained. In the light-colored region: self- 
sustained oscillations. Initial conditions: {Sx)t=o = 1/2 
and {Sy)t^o = {Sz}t^o = 0, {Jx}t=o = {Jy}t^o = j/2 and 
{Jx)t=o = j/V2. 



tions can be obtained through the following expression: 



(n.(i)) = i((iV(0))e-" + ^^^l±lH(i 



(7) 



which relates the level occupation with the z-component 
of the electron spin. If the coupling between the elec- 
tron and the large spins is isotropic (A^^ = Aj^ = A^), it 
is straightforward to see that in the stationary limit the 
spins decouple, and the well known Fano- Anderson solu- 
tion is obtained (see Appendix C). In contrast, we show 
below that the situation is drastically different for the 
anisotropic case where the stationary solutions for the 
EOMs depend on the coupling between the spins. 

The average electron current (/) through the QD is 
solely due to a decay at rate Tr^ from the spin- up QD 
level into the right lead. 



{i{t)) = eFHt(nt(t)), 



(8) 



where e denotes the electron charge. In the long-time 
limit of the current can be written as (see Eq. ([7])) 



eF 
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(9) 



Henceforth, for convenience we take F^f = Tr^ = F/2. 
Other options give similar behaviour except for the tran- 
sient solutions. 



III. REGIONS IN PARAMETER SPACE 

The stationary solutions of the EOMs, Eq. ([6]), can be 
obtained analytically and we find eight fixed points. Two 
of these fixed points, however, always have a finite imag- 
inary component, and as they have no physical meaning. 



we leave them out of the subsequent analysis. The re- 
maining fixed points serve to divide the parameter space 
of the model into distinct regions, as shown in Fig. [2j 
Introducing the notation 

the six relevant fixed points are 

P± = (^0, 0,-^0,0, , (11a) 

= (^0,62,-:^,|-S2,±Si,-^) , (lib) 
PlL2± = (^0,-B2,-^,-^B2,±B^,^^^ ,(llc) 
where 




For certain values oi B^, F and A, the quantities Bi and 
B2 f (|12ap and Eq. (jl2bp ) can have finite imaginary com- 
ponents and therefore, points 7^ii,i± and 7^ii,2± only have 
physical meaning in the region of parameter space where 
Bi and B2 are real. Fig. shows a projection of the 
3-dimensional parameter space on the F versus plane 
for a fixed A. This diagram is divided in three regions. 
In region I, Bi is a pure imaginary number, and hence, 
Vii^ij^ and 7'ii,2± are nonphysical, and V± the only phys- 
ical fixed points. In region II, Bi and B2 are both real, 
and all six fixed points are physical. In region III, B2 
is purely imaginary, and again V± are the only physical 
fixed points. Points 'P± are thus physical solutions for 
the EOMs Eq. ([6]) in all three regions, whereas the fixed 
points Pii,± are physical only in region II. The boundaries 
between the regions are obtained by solving the equations 
Bi —Q and B2 — Q, namely 



61 = ^ F : 



i f - {Bjxy B, 

l/A-B,J\ A 



^2 = 0^^ = 1, (13) 
and these two equations give the lines plotted in Fig. ^Bji. 

A. Region I 

In order to obtain the time evolution of the electron 
and large spin components and the electronic current 
through the QD, the EOMs Eq. ([6]) are solved numeri- 
cally. Fig.|3]shows the time evolution of the electron spin 
and the large spin components, and the current through 
the QD in region I of the parameter space. All exhibit 
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FIG. 3. (Colour online.) Time evolution in region I of a) 
the electron spin components, b) the large spin components, 
and c) the current through the QD obtained by solving nu- 
merically the EOMs of Eq. ((6]). In this region, the solutions 
exhibit a slow damped behaviour. In the stationary limit, the 
large spin is completely polarized in the direction parallel to 
the external magnetic field, and a QD electron trapped in the 
spin-down state. Current is due only to tunnelling through 
the spin up level and in the stationary limit tends to a con- 
stant value of 5/8. The parameters here are B^/X — 0.1 and 
F/A = 9, with initial conditions {Sx)t=o = 1/2, {Sy}t=o = 
{S,)t=o = 0, {X)t =o = { Jy)t=o = (5/y2)(y5 - l)/2 and 



completely damped oscillations. In the previous discus- 
sion, we have seen that in region I, 'P± Eq. (Illap are the 
only physical fixed points. Depending on the choice of pa- 
rameters and initial conditions, the system will evolve to 
V+ or V- ■ For the parameters and initial conditions cho- 
sen in Fig. [3l the system evolves towards the fixed point 
V+. In this case, the large spin becomes completely po- 
larized in the direction parallel to the external magnetic 
field (Fig.[5]D), and a spin-down electron remains trapped 
in the QD (Fig. [3^) and the interaction between the elec- 
tron spin and the large spin is no longer effective. Spin- 
up electrons, however, can still tunnel through the QD 
(Fig.|3t), and in the stationary limit the current becomes 
(see Eq. ^) 



iHt)) 



eF 



(14) 



In region I, then, the coupling of the two spin systems 
with the external leads results in complete damping of 
the transient oscillations of the electron and the large 
spin components and the current. A finite, fully spin- 
polarized electron current flows through the QD that in 
the stationary limit is not influenced by the interaction 
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with the large spin. 



B. Region II 



In region II, the EOMs Eq. © exhibit both damped 
and self-sustained oscillatory solutions, depending on the 
choice of parameters and initial conditions. Fig.[2}3 shows 
the part of region II where the self-sustained oscillations 
are found. This behaviour can be seen for all intensi- 
ties of the external magnetic field in region II, but only 
for small values of coupling F with the leads. Compar- 
ing Fig. [2^ and Fig. [Jb we see that most values for Bz 
and F in region II lead to damped oscillations. Further- 
more, although we have given analytical expressions for 
the boundaries between the different regions Eq. ([T3| . 
we have not found an expression for the boundary be- 
tween the regions inside region II where self-sustained 
and damped oscillations are found. Fig. has been ob- 
tained by solving the EOMs Eq. ([6|) in region II. As can 
be seen, the boundary between both regions is fuzzy in 
contrast with the ones obtained between regions I, II and 
III Eq. ([Ta]). Moreover, Fig.[2lD shows small "islands" in 
the oscillatory region, where damped solutions are ob- 
tained. 



1. Damped Oscillations 

Fig.H^, Fig. lib and Fig.Hl;, show the time evolution of 
the electron and the large spin components, and the cur- 
rent in region II with parameters Bz and F such that they 
all exhibit damped oscillations. Previously we have seen 
that in region II all the six fixed points are physical. For 
the parameters and initial conditions chosen in Fig. |4^, 
Fig. and Fig. HJ;, the system evolves towards the fixed 
point Vii.i+. The large spin becomes almost completely 
polarized in the y-direction (Fig. |41d), perpendicular to 
the external magnetic field, and the current becomes 
(see Eq. (El)) 



eF 



3 

4 



X 



(15) 



Thus, the stationary current increases if either the exter- 
nal magnetic field decreases or the coupling between the 
spins increases. Since in region II Bz/X < 1/4, the cou- 
pling between the electron and the large spins enhances 
the current through the QD, compared with the current 
obtained in region I (Eq. ([Ti)) ). Nevertheless, the result of 
coupling the two spins to the leads stills yields complete 
damping of both spin oscillations, as in region I. 



2. Self-Sustained Oscillations and Chaos 

We shall now focus on the small region in region 
II where self-sustained oscillatory solutions are found 



(Fig. [Id). Fig. Fig. gj? and Fig. gf, show the time 
evolution of the electron and the large spin components, 
and the current through the QD. The chosen values of 
Bz and F lead to complicated, but periodic, undamped 
oscillations. Fig. [5^ shows the Fourier spectrum of the 
current time evolution of Fig. gf in the long-time limit. 
The spectrum exhibits peaks at well defined frequencies, 
which clearly confirms the periodic behaviour of the cur- 
rent. Furthermore, in non-linear systems, self-sustained 
oscillations are a signature of limit-cycles and in Fig. [5^ 
and Fig. [5}d we plot the electron and the large spin tra- 
jectories in phase space, projected on the x-z plane, in 
the long-time limit. These figures show that the spin 
trajectories are precisely limit-cycles. For all the ini- 
tial conditions chosen, the system always converges to 
them. Finally, Figs. |3^, |3)i and |3J show that decreasing 
F turns the periodic self-sustained oscillations chaotic. In 
this case, the Fourier spectrum of the current, shown in 
Fig. [5}d, is uniformly distributed through all frequencies, 
which is a clear signature of chaos. Fig. [St and Fig. [SJi 
show the electron and large spin trajectories in the long- 
time limit, where it can be seen that they perform com- 
plicated non-periodic paths. In this area of region II, 
the coupling between the interacting spins and the leads 
does not produce damping of the spins as in the previous 
cases. Moreover, the electron current through the QD 
captures the complicated dynamics due to the interac- 
tion between the electron and the large spin, as seen in 
Fig. \S and Fig. [1. 



C. Region III 

Fig. [5^, Fig. [5)d and Fig. [5J; show the time evolution of 
the spin components and current for typical parameters 
in region III. They all exhibit periodic self-sustained oscil- 
lations. Fig. [5] shows the different limit-cycles performed 
by the electron spin in phase space, projected in the y-z 
plain, when the value of the external magnetic field is in- 
creased. The trajectories found for the large spin in the 
long-time limit suggest that this behaviour can be under- 
stood by means of an effective model in which the large 
spin simply acts on the QD electrons as an ac magnetic 
field in x-direction with amplitude 

Bac(i) = ^ icos{Bzt) - sm{Bzt)) . (16) 

The EOMs for this effective model are (see Appendix D 
for details) 



~Bz{Sy) 

Bz{Sx) - 



Bac{t){Sz 



(17) 



Thus, in this region the six autonomous non-linear equa- 
tions Eq. ^ can be approximated by a set of three non- 
autonomous linear equations Eq. (1171) . The agreement 
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FIG. 4. Results for three different parameter sets are shown, each of which gives rise to very different system behaviour. Figures 
a), b), c) show fast damping behaviour (-Bz/A = 0.2, F/A — 0.7, dark region in Fig. [2}5). In the long-time limit the large spin is 
almost completely polarized in the direction perpendicular to the external field, but unlike in region I, the spin down electron can 
escape from the QD into the right lead due to the interaction with the large spin. Figures d), e), f) show periodic self-sustained 
oscillations (Bz/A = 0.1, F/A = 0.16, light region in Fig. ^jp), which is a signature of limit-cycles in phase space (see Fig. [S]). 
Figures g), h), i) show chaotic self-sustained oscillations {Bz/X = 0.1, F/A = 0.015, light region in Fig. [2]d). In the oscillatory 
cases, the oscillations are captured in the current through the QD, and in particular, the chaotic behaviour is observed in the 
current (Fig.H). T he initi al conditions are {S^)t=o = 1/2, {Sy)t=o = {Sz)t=o = 0, {Jx)t=o = {Jy)t=o = (5/V2)(V5 - l)/2 and 
{Jz)t=o = (5/^/2)^/5 + ^/5. 



between the solutions obtained with this effective model 
and the full EOMs is very good. 

In region III, the coupling between the two spin system 
leads to self-sustained oscillations which are visible in the 
electron current through the QD, as shown in Fig. Ht. 



IV. CONCLUSIONS 

We have studied electron transport through a quan- 
tum dot coupled to ferromagnetic leads, in which the 
electron spin interacts with a large spin while an exter- 
nal magnetic field is applied. We have found that the 
motion of the electron spin, the large spin and the cur- 
rent through the QD strongly depend on the coupling be- 
tween spins. When the electron spin and the large spin 
are isotropically coupled, the large spin becomes com- 
pletely polarized and decouples from the electron spin. 
Conversely, when the electron spin and the large spin are 
anisotropicaly coupled, we have found that their motion 
and the current through the QD can either behave as 
in the isotropic case or show self-sustained oscillations 



which, furthermore, can be periodic or chaotic. Switch- 
ing between different behaviours can be obtained by vary- 
ing either the strength coupling with the leads or the 
intensity of the external magnetic field. 

We foresee two possible experimental realisations of 
the large spin of our model. The first is as an effective 
model on a hyperfine bath. Here a semi-classical treat- 
ment may be justified by considering that the number of 
nuclei spins in semiconductor QDs interacting with an 
electron spin is very large (e.g., for GaAs QDs there are 
typically 10^-10^ nuclei spins). Situations in which the 
hyperfine interaction is anisotropic have been discussed 
in Refs. [ssl-fSTj. The second realisation is that our large 
spin represents the spin of a magnetic impurity of a doped 
semiconductor or a magnetic atom in a single molecular 
magnet. While in this case the spin may not be so large, 
mean-field analyses such as pursued here can still provide 
useful information, e.g. Ref. [ssj . 

From the theoretical point of view, it would be inter- 
esting to investigate how the features of this semiclassical 
treatment are reflected in a quantum master equation ap- 
proach, in which the electron and the large spin are both 
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FIG. 5. Fourier spectra of the non-damped current time evo- 
lutions shown in Fig. 2] in the long-time limit. Figs, a) and 
b) show the Fourier transform of Fig [If and|4l, respectively, 
where v is the frequency. Figure a) shows peaks at well de- 
fined frequencies, meaning that behaviour of the current is 
periodic. However, figure b) shows a uniform frequency dis- 
tribution, which is a signature of chaotic dynamics. 
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FIG. 7. In this region the solutions exhibit periodic self- 
sustained oscillations, which are reflected in the current. The 
corresponding limit-cycles are shown in Fig. [7]). The pa- 
rameters chosen here are B^/X = 1.0, F/A — 10. The ini- 
tial conditions are {Sy)t=o = 1/2, {Sx)t=o = {Sz)t=o = 0, 
{Jx)t=o = {Jy}t^o = 3,/ll/2 and {Jz)t=o = -1. 
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FIG. 6. Electron spin (left figures) and large spin (right) 
trajectories projected on a two dimensional plane for the non- 
damped solutions in region II (light region in Fig. [2^). Figures 
a) and b) show the formation of a limit-cycle as seen in the 
time evolution plots. Fig. |4jl, Fig. |4^, and Fig. |4f (B^/X = 
0.1, F/A = 0.16). Figures c) and d) correspond to the time 
evolution plots of Fig. |4^, Fig. |4]i, and Fig. |4j which suggest 
that the trajectories are chaotic {B^/X — 0.1, F/A — 0.015). 




-0.2 0.2 0.4 -0.2-0.1 0.1 0.2 0.3 

2(4) 2(5,) 

FIG. 8. Electron spin trajectories projected in the {Sy)-{Sz) 
plane in region III. The solutions of the EOMs given in Eq. ([6]) 
are periodic self-sustained oscillations (Fig. [T]). Figures a)- 
f) show the different limit-cycles obtained when varying the 
external magnetic field. F/A = 1. 
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Appendix A: Derivation of the equations of motion 

In this appendix, we summarize the steps in the derivation of the EOMs (Eq. ©)• We start with the Hamihonian 
of Eq. (H]) and, for later convenience, we shift the reservoir frequencies J2ikcr <^ikac]j^^cika -> I]/fc<T (^ifco- + m)^\kcr^ika 
where /Uj is the chemical potential of lead /. Under the mean- field approximation considered in this work (see Eq. ([5])), 
the closed set of EOMs obtained for the time evolution of operators in the Hamiltonian (Eq. ^) are then computed 
to be: 

i ^ {did,,) - ^ (<5.'t(4rf\> + 5,'Md^) ~ 5,^{dld,') - 5ai{d\d,')) {X) 

+ i (^.-t (dUt) ~ ^-'l (dUi) - (5,t {d\da') + S^i (did,,)) (a(J,) + B,) 

- E {^i^i^lkM - likidlhk.')) (Al) 

I, k 

iJt) = \w{^ikA) + \ {mD + Br) {cij^) + 7r'fc'(cLci'fe't) - i^ika + w)(5Lrft) - itkidUt) 

l',k' 

■LJi) = ^WiclA) - I {MJz)+b) (clJi,) + it'k'icL^l'k'i) ~ ieika+^il){clJi)-ltk{dUl) 



l',k 

i 4l {c]i.„d_ " 
dt - 

l',k' 

and 



dt 

{Jy) = -X{S,){X) + (a(^,) + B,) (J,) 



dt 
d_ 
di 



{Jz)^X{Sa.){Jy) (A2) 



where we have used the choice = Xz — ^ and Xy = 0. Since the EOMs for the large spin components have already 
the desired form (see Eq. ©), hereinafter we shall focus on the time evolution of the electron operators (Eq. (|Aip ). 
Under the Born approximation the leads are assumed to be in thermal equilibrium for all time, 

{^Ika^l'k'o') = fla Sw Saa' 6{k' - fc), (A3) 

with fi„ the equilibrium Fermi-Dirac distribution for spin-cr electrons in lead /: 

= /(,,,,) = ^^^^^J^^^^^ - (A4) 
Applying the Laplace transform, {A)s = /q°° e'^'^*{A)t dt, to Eq. (jAip we obtain: 

is {dlda')s = {5a'^{d\di)s + 6a'i{dld\)s - 5a^{d\da')s - 6ai{d\da')s^ {Jx)s 

+ \ (-^^'t {dld^)s - Sa'i {didi,)s ~ (5<xt {d\da')s + {dld,')s) (a(J;,), + B^) 

- E (7/fe(cL<,4')s - 7rfc(4c/fc<T')s) + i {dU<r')o (A5a) 



and 



{cjkadt)s = ^{Jx)Acjkadi)s + ^ (A(J^)s + b) (4fca^t)« + /'^ '^^t Tifc " i<^lka + tM){c\kA) s " llk{dld^) k^^^) 

){Jx)s{c\kA)s - ^ (a(J^), + b) {c\f,Ji)s + fla 5ai 7ik - {<^ika + IJ'i){c\kJi) s " 7zfe (4^^^) jJA5c) 
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where {A)o denotes the expectation vahie of operator A at time t = 0, and where we have taken {clj^^dcr')o = 0. After 
some algebra, Eq. (jA5bp and Eq. (|A5cl) become: 



{fio- Sat - (dldt) 



2j^,X{X)s{{dUl)s~.fi.) 



{!ia5at-{did^),) 



(A6a) 



2itkHJ.)s{{dldt)s-fia) 



(A6b) 



We now consider the infinite bias hmit and set, for the left lead, ^ oo, and for the right, /i/j — > — oo. In this limit, 
the denominator of the first term in Eq. (jA6al) becomes eika + /i/ + iO+, with positive infinitesimal 0^, and the second 
term is seen to be of the order ^^"^ and thus negligiable compared with the first term (of order j-i^^). Equations (IA6a[) 
and (|A6b[) thus become: 



e;feo- + Mi + 



(A7) 



This result allows us to rewrite the summation that appears in Eq. (|A5a[) as: 



Ik I 



e + M/ + «0+ e + fj.1- iO+ 



{f{e)Saa' - {dlda,)s) 



with the lead- and spin-dependent rates 



(A8) 



with piaif) density of states of the /-th lead. We assume these rates to be energy-independent, r;cr(e) = ^la (flat band 
approximation). Using the Sokhatsky-Weierstrass theorem, 

=p i7r(5(x). 



upon evaluation of the Fermi functions at /i^ = oo and pn = —oo, we obtain 



E 

Lk 



llk{4kad'a')s - l!k{dlcika')s] = ;j E^^''^ + ^ la') {dlda') s " i^bJa 



(A9) 



Replacing the previous expression in Eq. (jA5al) gives: 

{hcr)s = X{J^)s{Sy)s {Sat - ^ai) - T^{na)s + ^ La 



S {Sy)s = -X{S,)s{X)s + (5,), [X{X)s +B,j-r {Sy)s + {Sy)o 

s {S,)s = X{.l),{Sy)s - T{S,)s + i (Fit - Tlx) 
where Ta —TLa + ^Ra for a =ti4-: although we have assumed for simplicity r-|- = Fj, = F; and the identities: 



(AlO) 



Sx — 


\ {d\di - 




Sy = 




-did,) 









(All) 
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have been used. Finally, inverse Laplace transforming Eqs. (jAlOp yields the EOMs ([6]) for the occupation and the 
spin components of the electron in the QD. 

The EOM for the total occupancy of the QD is obtained by summing the EOMs of the spin-up and spin-down 
occupations, 

±{N)=-T{N)+rL^ + TLi. (A12) 
Notice that this EOM is independent of the electron and large spins, moreover, it is exactly solvable giving; 

(Nit)) = (iV(0)) e-" + ^ ^^-^ (1 - e-") . (A13) 



1. Heuristic derivation 



The electronic part of our EOMs can be seen to make sense by considering a more intuitive derivation using rate 
equations for the QD occupations (Eq. ([6])) when A = 0. The QD states in the transport window are {|0), | f): \ i), 
|2) = I ti)}- In the infinite bias regime electrons can tunnel into the QD from the left lead and tunnel out of the QD 
to the right lead, thus 

Po = M^o>t + W„1pi - (Wli, + Wl'„)po 

h = KoPO + W^2P2 - {W^i + W^{)p^ 

P2 = M^sVt + W^2>i - ( W^t'2 + )P2 . ( A14) 

where pi is the probability of finding an electron in state \i). Wj^ is the tunneling rate from the initial state |i) to the 
final state |/) through the ^-th barrier. Using the conservation of total probability (Tr(p) — 1) we get 

= W^t^o(l - Pi) + (Wf2 - W^t^o)P2 - (W^t^o + W^ot + Wi\)p^ 

Pi = - pt) + (wf2 - wio)P2 - (w^i, + + wi^jp^ 

P2 - M^2tPt + W^2>i - (W^t'2 + W^2)P2- (A15) 

We now consider that W^^^ = W^l = and W^^^ = W^^ = Fl^, and W^^ = W^^ = F^t and Wg^ = W^^ = F^^^, so 

Pt = rLt(l - Pi) + i^Rl - ^Lt)P2 - (FLt + ^Bt + ^Li)Pt 
Pi = rL|(l - Pt) + {^m - ^Li)P2 - (FLi + ^Ri + ^Lt)Pl 

P2 = ^LlPt + TitPi - (r^t + ^Tfi)P2- (A16) 

Finally, since 

Pa +P2= Flo- - (Flo- + F_Ro-)(pc -I- P2), (A17) 

and (ha-) — Pa + P2, we arrive to 

(A,) =-F(fv)+Fi, (A18) 



where we have used that F^o- + ^Ra — F. 
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Appendix B: Isotropic model 

The EOMs for the completely isotropic 

A (n.) = A [{l){Sy) - {Jy){S,)) (<5,t - - nna) + T 
I (S,) = X{Jy){S..) - [X{1) + B,) (Sy) - r{S,) 
I (Sy) - -X{X){S,) + (a(J,) + B,) (5,) - r{Sy) 
^ (5.) = A [{l){Sy) - {Jy){S.)) - r{S.) + i (Fit - Tl;) 
^ (J,) = A(4)(J,) - (a(5,) + S,) (J,) 
I (J,) = -A(5,)(J,) + (a(^,) + S.) (J,) 

I ( J.) = A ((5.) (J,) -(5,) (J.)). (Bl) 

To find the solutions in the stationary hmit we put to zero the time derivatives. Therefore, it can be seen right away 
that in the long time limit the quantum dot occupations decouple from the large spin components and becomes 



(B2) 



Thus, the spin dynamics can not be observed in the current. 



Appendix C: Effective model for region III 



In this appendix, we summarize the steps in the derivation of the effective EOMs (Eq. (|17p ) for region III of the 
parameter space (Fig. [2^). Applying the transformation (S) = e^^* TZ{t) ■ (S) and (J) = TZ{t) ■ (J) with 



cos{B,J) -sm{B^t) 
TZ{t) = \ sm{B^t) cos(B^t) 
1 



(CI) 



to the EOMs (Eq. ([6])) they become: 



^ (S.) = -x[{J,){Sy) + [(X) cos{B,t) - (Jy) sm{B,t)\ (^,) sm{BJ)} 
^ (Sy) = X {(J.) (5,) - [(J,) cos{B,t) ~ {Jy) sin(B,i)] (5,) cos(B,t)} 



dt 
d 

It 
d 

di 
d_ 

dt 
d 

di 
d_ 

dt 



(S,) = A {J^){Sy) cos^B.t) - {Jy){S,) sin^B^t) + i{J^){S^) - {Jy){Sy)) sm{B,t) cos{B,t) + -(E^t - J e 



(Jcc) = -Xe-" [{S,){Jy) + [{Sec) cos{B,t) - (Sy) sin(B,t)] (J,) sm{B,t)} 
(Jy) = Xe-^' {{S,){X) - [{S^) cos{B,t) - {Sy) sin(S,i)] (J,) cos(B,t)} 

{S,){Jy)cos'{B,t) - {Sy){J.,) sin^B.t) + ((^,)(J,) - {Sy){Jy)) sm{B ,t) cos{B ,t) 



J.) = Xe 



(C2) 



Since in the long-time limit d (Ji) /dt — > 0, we assume (J) to be stationary. Therefore, the EOMs for the electron spin 
in the original frame become: 



B,tf X (s) - r(s) 



1 



(Fit - Tl;) u. 



dt' ' ' ^ ' ■ 2 ' 

where Boff = (-Bac(0: 0, B^) and is the unit vector pointing in the z-direction. 



(C3) 
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